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Tactiele  displays  en  elastische  golven 


Tactiele  displays  zijn  displays  voor  de  huid,  zoals  beeldschermen  visuele  displays  zijn  voor  de 
ogen.  Een  bekend  voorbeeld  is  de  trillende  mobiele  telefoon.  In  de  mobiele  telefoon  zit  een 
trilmotortje  dat  aangaat  bij  een  binnenkomend  gesprek.  TNO  Defensie  en  Veiligheid  heeft  een 
tactiel  display  in  een  vest  met  wel  64  van  deze  trilmotortjes,  ook  wel  tactors  genaamd. 


Probleemstelling 

Om  op  zo’n  groot  display  informatie 
duidelijk  af  te  beelden  moet  ervoor  gezorgd 
worden  dat  de  trillingen  van  de  tactors  zo 
lokaal  mogelijk  op  de  huid  worden 
aangeboden.  De  trillingen,  ook  wel 


elastische  golven  genoemd,  mogen  zich  met 
te  veel  verspreiden  door  de  huid  of  door  het 
materiaal  waarin  de  tactors  zitten  bevestigd. 
In  het  huidige  tactiele  vest  gebeurt  dit  wel. 
In  opdracht  van  het  Ministerie  van  Defensie 
is  een  Nationaal  Techno logie  Programma 


gestart,  waarin  demonstrators  zullen 
worden  gemaakt  voor  een  nieuwe  generatie 
tactiele  displays.  Een  onderdeel  hiervan  is 
het  opstellen  van  ontwerpcriteria  voor 
tactiele  displays,  die  de  verpreiding  van 
trillingen  minimaliseren. 

Beschrijving  van  de 
werkzaamheden 

De  ontwerpcriteria  zijn  opgesteld  met 
behulp  van  een  literatuurstudie  waarin  de 
generieke  theorie  van  elastische 
golfgeleiding  in  materialen  is 
gecombineerd  met  kennis  over  de  structuur 
van  de  huid. 

Toepasbaarheid 

De  ontwerpcriteria  die  uit  deze 
literatuurstudie  volgen,  zullen  worden 
meegenomen  in  het  ontwerp  van  een 
nieuwe  generatie  tactiele  displays  en  ervoor 
zorgen  dat  informatie  beter  kan  worden 
overgebracht  op  de  huid. 
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Samenvatting 


Probleemstelling 

Wanneer  tactoren  van  een  tactiel  display  in  een  drager,  zoals  een  vest  of  een  stoel,  zijn 
gemonteerd  spreiden  die  vibraties  zich  uit  door  de  huid  en  de  omgeving  wat  tot  gevolg 
heeft  dat  de  trillingen  niet  meer  lokaal  worden  waargenomen  en  moeilijker  zijn  te 
onderscheiden.  Vibraties  in  een  dergelijk  systeem  moeten  lokaal  zijn  voor  een  goede 
detectie;  spreiding  door  de  huid  of  de  omgeving  moet  geminimaliseerd  worden. 

Werkwijze 

Er  is  een  literatuurstudie  gedaan  die  een  kwalitatieve  beschrijving  van  de  tactor-skin 
interface  geeft.  Generieke  golfgeleidingstheorie  en  kennis  van  de  structuur  van  de  huid 
zijn  gecombineerd  om  ontwerperiteria  voor  tactiele  displays  te  bedenken  die  de 
verspreiding  van  trillingen  minimaliseren. 

Resultaten 

Volgens  de  literatuur  kunnen  drie  soorten  golven  zich  voortplanten  in  de  huid  en  de 
drager  waarin  de  tactoren  zitten:  schuifgolven,  drukgolven  en  oppervlaktegolven.  Zij 
planten  zich  radieel  voort  met  verschillende  snelheden  en  met  verschillende  golflengtes. 
De  golven  spreiden  zich  en  doven  uit,  maar  ze  kunnen  ook  energie  verliezen  door 
frictie.  Elke  type  golf  dooft  met  een  andere  snelheid  uit. 

Conelusies  en  aanbevelingen 

Tactoren  die  een  minimale  spreiding  van  trillingen  door  de  huid  of  de  omgeving  hebben 
trillen  parallel  aan  het  huidoppervlak  met  een  frequentie  van  200  Hz.  De  tactoren 
zouden  aan  de  huid  gelijmd  moeten  zijn  of  gecoat  met  een  materiaal,  zodat  de  tactor 
niet  over  de  huid  kan  glijden.  Het  zou  nog  beter  zijn  als  er  een  materiaal  met  een  dikte 
van  een  kwart  golflengte  tussen  de  tactor  en  de  huid  wordt  geplaatst.  Dit  zorgt  voor  een 
goede  overdracht  van  vibratie-energie.  Wanneer  de  tactoren  niet  slipvrij  op  de  huid 
bevestigd  kunnen  worden,  dan  is  het  beter  verticaal  op  de  huid  te  trillen.  De  grootte  van 
het  contactoppervlak  van  de  tactor  met  de  huid  moet  klein  zijn.  Tussen  de  tactor  en  de 
drager  zou  een  zacht  materiaal  met  daaromheen  een  stijf  materiaal  geplaatst  moeten 
worden,  zodat  er  zo  weinig  mogelijk  trillingsenergie  wordt  overgebracht  op  de  drager 
en  zoveel  mogelijk  energie  naar  de  huid  wordt  gereflecteerd.  Als  laatste  zou  de  grootte 
van  trillende  oppervlakten  van  bijvoorbeeld  de  behuizing  van  de  tactor  minimaal 
moeten  zijn,  om  ervoor  te  zorgen  dat  de  tactoren  geen  lawaai  maken. 
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Summary 


Purpose 

When  tactors  of  a  tactile  display  are  mounted  in  a  tactor  carrier,  for  example  a  vest  or  a 
chair,  vibrations  will  spread  out  through  the  skin  and  the  tactor  carrier.  As  a  result  the 
localisation  of  the  vibration  sources  is  reduced.  Vibrations  in  such  a  system  should  be 
local  to  ensure  good  detection;  spreading  through  the  skin  or  the  environment  has  to  be 
minimised. 

Method 

A  literature  review  that  gives  a  qualitative  description  of  the  tactor-skin  interface  has 
been  carried  out.  Genera!  wave  propagation  theory  and  knowledge  about  the  structure 
of  the  skin  haven  been  combined  in  order  to  create  several  design  criteria  that  minimise 
the  spreading  of  vibrations. 

Results 

From  the  literature  it  became  clear  that  three  types  of  waves  can  propagate  through  the 
skin  and  the  tactor  carrier:  shear  waves,  compression  waves  and  surface  waves.  They 
propagate  through  material  in  a  spherical  fashion  at  different  propagation  velocities  and 
with  different  wavelengths.  The  waves  spread  out  and  dampen,  but  they  can  also  loose 
energy  by  friction.  Each  type  of  wave  dampens  at  a  different  rate. 

Conclusions  and  Recommendations 

Tactors  that  have  a  minimal  spreading  of  vibrations  through  the  skin  or  the  tactor  carrier 
oscillate  parallel  to  the  skin  and  use  a  stimulation  frequency  of  200  Hz.  The  tactors 
should  be  glued  to  the  skin  or  coated  with  a  material  that  does  not  slip  easily.  It  is  even 
better  to  insert  a  material  with  a  thickness  of  a  quarter  of  a  wavelength  between  the 
tactor  and  the  skin.  This  results  in  a  better  transfer  of  vibration  energy.  When  the  tactors 
cannot  be  attached  free  of  slip  to  the  skin,  it  is  better  to  use  tactors  that  oscillate 
perpendicular  to  the  skin.  The  size  of  the  contact  area  between  the  tactor  and  the  skin 
should  be  small.  Soft  material  encased  in  stiff  material  should  be  used  between  the 
tactor  carrier  and  the  tactor  in  order  to  reduce  the  amount  of  vibration  energy  emitted  to 
the  tactor  carrier  and  to  reflect  energy  to  the  skin.  Finally,  the  area  of  vibrating  surfaces 
of  e.g.  the  housing  of  the  tactor  should  be  kept  minimal  to  reduce  the  noise  production 
of  vibrating  tactors. 
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1  Introduction 


The  past  six  year  TNO  has  investigated  how  information  can  be  displayed  on  the  skin. 

It  turned  out  that  an  array  of  vibrating  elements,  so-called  tactors,  placed  close  to  the 
skin  can  be  used  as  a  tactile  display.  TNO  has  built  several  versions  of  a  vest  with 
64  tactors  and  a  belt  with  eight  tactors.  They  were  used  in  several  experiments  that 
showed  that  it  was  possible  to  provide  the  wearer  of  a  tactile  display  with  enough 
information  to  enable  him  or  her  to  successfully  hover  with  a  helicopter,  to  navigate 
along  way-points  and  to  recover  from  spatial  disorientation  [Van  Erp,  Veltman,  Van 
Veen,  &  Oving,  2002;  Van  Erp  &  Van  Veen,  2003;  Van  Erp,  Jansen,  Dobbins,  &  Van 
Veen,  2004;  Van  Veen  &  Van  Erp,  2003;  Van  Veen,  Spape,  &  Van  Erp,  2004], 

The  current  hardware  of  the  tactile  display  is  designed  to  work  in  laboratorial 
circumstances  and  it  is  not  very  suitable  for  field  experiments.  In  the  project  called 
Tactile  Displays  (TACT)  new,  state-of-the-art  hardware  that  can  be  integrated  in 
garment  is  developed.  It  is  designed  to  be  used  in  the  field  and  still  provide  a  clear 
tactile  sensation.  Several  demonstrators  are  designed  that  can  show  the  capabilities  of 
tactile  displays  on  land,  in  the  air  and  in  the  water. 

In  order  to  provide  clear,  local  vibratory  stimulation  to  the  skin,  the  spreading  of 
vibrations  through  the  skin  or  the  garment  has  to  be  minimised.  This  report  is  a 
literature  review  that  gives  a  qualitative  description  of  the  tactor-skin  interface  when  the 
tactors  are  mounted  in  for  example  a  tactile  vest.  It  combines  knowledge  of  wave 
propagation  in  materials  given  in  Chapter  2  with  knowledge  of  the  structure  of  the  skin 
and  of  its  mechanicals  properties  given  in  Chapter  3.  In  Chapter  4  design  criteria  for  a 
general  tactile  display  are  composed. 
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2  Theory  of  Wave  Propagation  in  Materials 


2.1  Mechanical  Properties  of  Materials 


The  propagation  of  waves  in  a  materia)  is  strongly  related  to  the  mechanical  properties 
of  that  material.  Two  important  properties  are  the  elasticity  and  the  density  ( p )  of  the 

material.  When  a  uniform  uniaxial  stress  (force  per  surface  area) 

[N/m2]  (1) 


is  applied  to  a  rod,  a  change  A  /  in  the  length  /  of  the  material  will  occur.  The  relation, 
according  to  Hooke’s  law,  is 


e 


X 


A/ 

1 


=  ^ -  [N/m2], 


(2) 


in  which  E  is  called  the  Young’s  modulus  and  £  the  strain.  It  is  obvious  that  when  the 
Young’s  modulus  is  high,  the  material  hardly  changes  in  length  when  stress  is  applied. 
The  material  thus  has  a  high  stiffness.  Most  materials  and  certainly  the  skin  show  a 
nonlinear  relationship  between  the  stress  and  the  strain.  The  parameters  are  often 
linearised  for  small  changes  in  a  or  £  .  Not  only  does  a  change  in  stress  cause  a 
change  in  the  length  of  the  material,  but  also  a  change  in  the  width  of  the  material  (both 
in  y  and  z  direction);  in  the  case  of  an  elongation  of  the  material,  the  width  will 
decrease.  The  ratio  between  the  relative  change  in  width  and  length  is  called  the  Poisson 
ratio  and  is  defined  in  an  isotropic  medium  as 


Ad/d  _  8 y  _  £, 

A ///  £ x  £x  ’ 


-1  <  v  <  — 

2 


(3) 


in  which  Ad/d  is  the  relative  change  in  width.  In  most  materials,  v  is  between  zero  and 
a  half.  For  v  =  1/2  the  material  is  incompressible  and  for  v  =  0  the  material  is  fully 
compressible.  When  Hooke’s  law  is  applied  to  a  cube,  the  strain-stress  relation  in  three 
dimensions  becomes: 


V 

i 

1  -  V  -  V 

£„ 

1 

-v  1  —  V 

y 

£z 

—  V  —v  1 

The  inverse  of  Equation  4  is  given  by 
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2.2 


(7  = 


(l  +  vXl-2v) 


1  —  V 

V 

V 


X  = 


A  +  '2.fu  A 

A  A  +  2/j 
k  A 

vE 


V  V 

I  -v  V 

V  1  -v 

/l 


■e 


A 

A  +  2jj 


■s,  with 


(5) 


in  which  A  and  ju  are  called  the  Lame  constants.  These  two  are  often  used  in  the 
literature  to  describe  the  elastic  properties  of  a  material,  ju  or  G  is  also  known  as  the 
rigidity  of  a  material,  the  shear  modulus,  or  the  shear  elasticity;  when  /u  is  large  the 
material  will  hardly  change  in  width  when  uniaxial  stress  is  applied.  A  is  also  known  as 
volume  compressibility.  Basically,  A  indicates  how  easily  particles  in  the  material  can 
be  pulled  apart  or  pressed  together,  whereas  fd  indicates  how  easily  these  particles  can 

slide  along  each  other  [Weisstein,  2004b;  Leijendeckers,  Fortuin,  van  Herwijnen,  & 
Leegwater,  1998]. 


Longitudinal  Waves  and  Shear  Waves 

The  mechanical  properties  mentioned  above  are  important  for  understanding  the 
propagation  of  two  main  types  of  elastic  waves  that  can  travel  through  materials: 
longitudinal  waves,  also  known  as  P-waves  or  pressure  waves  or  body  waves  and  shear 
waves,  also  known  as  S-waves  or  transverse  waves.  In  the  first  case,  the  oscillation  of 
particles  in  the  material  is  in  the  direction  of  the  wave  propagation  (see  Figure  1).  In  de 
latter  case,  the  oscillation  is  perpendicular  to  the  direction  of  the  wave  propagation  (see 
Figure  2). 


wave  propagation  direction 

< - 


* - ► 

wavelength 

Figure  1  Schematic  presentation  of  a  plane  longitudinal  wave  that  propagates  through  a  material.  The 
particles  in  the  material  move  back  and  forth. 
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wave  propagation  direction 


wavelength 

Figure  2  Schematic  presentation  of  a  plane  shear  wave  that  propagates  through  a  material.  The  particles 

in  the  material  move  up  and  down. 

The  wave  equation  for  plane  longitudinal  waves  propagating  in  the  x-direction  is  given 
by 


^  =  ci 


. .  E  1  -  v  A  +  2ij 

,  Withe,  = -7 - T7 - — r  =  — -  , 

V  P  (l+vXl-2v)  V  p 


in  which  ux  is  the  displacement  of  particles  in  the  material  in  the  direction  of  the  x- 
axis,  c,  the  propagation  or  phase  velocity,  and  p  the  density  of  the  material.  The 
propagation  velocity  increases  monotonically  for  v  >  0 .  The  wave  equation  for  plane 
shear  waves  is  the  same,  except  for  the  propagation  velocity  cs : 


^:  =  c,>^,withc,=  E=  £ 

dt  dx  V  P  V  P  2(l  +  v) 


The  displacement  of  the  particles  is  now  in  the  direction  of  they-axis.  The  velocity  of 
shear  waves  is  always  smaller  than  the  velocity  of  body  waves  and  decreases  with  v  for 
v  >  0 .  Shear  waves  cannot  propagate  through  fluids  and  gasses,  because  there  is  no 
shear  coupling  between  particles.  The  general  solution  of  the  wave  equation  for  a  wave 
created  by  sinusoidal  excitation  and  propagating  in  the  positive  x-direction  is 


=  c 2  =>  u(t,x)  =  A -e^6*  kx\  with  k  =  —  =  —  and  a>  =  2nf . 

dt2  dx2  V  1  Ac 


A  is  the  amplitude  of  the  wave,/is  the  wave  frequency,  c  is  the  propagation  speed,  k  is 
called  the  wave  number,  and  A  =  c/  f  is  in  this  case  the  wavelength  [Weisstein,  2004c; 
Main,  1993;  Pain,  1993]. 

When  the  wave  propagates  through  the  material  a  dynamic  force  acts  on  each  particle  in 
that  material.  When  on  a  given  moment  the  force  on  a  particle  differs  from  its 
neighbouring  particle  there  will  be  local  stress.  The  difference  in  force  dF  acting  on 
two  particles  close  together  in  a  material  with  each  mass  m  is 

dF  da  d'u  .  ,  2  , 

—  ~m  —  -m — - — =  imkor  A-eA 

dx  dx  dt2dx  /rn 


=  mkco2  A-  e 


j(fOt-kx+7T  2) 


=  —  tv  A  ■  e 
c 


j((of-k x+z-2) 
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The  direction  of  the  force  is  along  the  axis  of  particle  motion.  The  maximum  amplitude 
of  the  force  difference  increases  with  the  wave  frequency.  Because  of  the  lower 
propagation  velocity  of  shear  waves  compared  to  compression  waves,  the  force 
difference  is  larger  for  shear  waves. 

2.3  Damping  of  Waves 


The  attenuation  of  elastic  waves  is  governed  by  two  mechanisms.  The  first  mechanism 
is  dispersion;  a  wave  driven  by  a  small  source  will  not  be  a  plane  wave,  but  a  spherical 
wave.  The  energy  of  a  spherical  wave  decreases  with  the  square  of  the  radius,  because 
the  area  of  the  sphere  increases  with  the  square  of  the  radius.  As  a  result  the  amplitude 
decreases  with  the  radius  [Feynman,  Leighton,  &  Sands,  1977;  Petrashen,  2003]: 

u{t,r)  =  -A-eA,0,-kr).  (10) 

r 


a  is  a  spherical  wave  propagating  outwards  and  r  is  the  radial  distance  from  the  source. 
The  factor  in  front  of  the  exponential  is  the  amplitude  in  which  A  is  constant. 

The  second  mechanism  is  conversion  of  wave  energy  into  heat  [Yerges,  1969].  This  can 
be  caused  by: 

-  the  viscosity  of  the  material; 

-  nonlinear  relations  (Hysteresis); 

-  frictional  damping  (Coulomb). 

Mathematically,  attenuation  is  introduced  by  making  the  Lame  constants  pi  and/or  A 
complex.  The  imaginary  parts  of  pi  and  A  are  called  shear  viscosity  and  volume 
viscosity  respectively.  As  can  be  seen  from  Equation  6  and  7  complex  Lame  constant 
results  in  a  complex  propagation  speed.  For  the  case  of  a  shear  wave  the  complex  wave 
propagation  velocity  c*  becomes 


.  •  .  A  pi  ,  fl 

P  =M  +JM  =>cs  =  —  =>  c,  =  — 
P  V  P 


1  + 


JJL 

2  pi 


when  pi  «  pi  .  (11) 


pi  is  often  made  linearly  dependent  of  the  wave  frequency.  Equation  1 1  leads  to  a 
complex  wave  number  k’ 


, .  co 

k  =  —  '■ 


co 


f 


r\ 
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\+JA 

v  2// 


ML 

2  pi 


=  k’~  jk 


(12) 


When  Equation  12  is  combined  with  Equation  8,  an  attenuation  term  appears  in  front  of 
the  solution  of  the  wave  equation,  which  is  dependent  of  the  position,  the  wave 
frequency,  the  density  of  the  material  and  the  rigidity  of  the  material: 


iy,x)  =  A-eJ^-k'x^  =  ^ 


-  A-e 


,4A 

in  T-n 


X  J\  (Ot-lO  I  — T.T 


(13) 


In  the  case  of  a  compression  wave  A  is  often  taken  real  and  pi  complex  in  the 
literature.  This  means  that  the  imaginary  part  of  k  is  smaller  than  it  would  be  in  the  case 
of  a  shear  wave  and  therefore  longitudinal  waves  experience  less  attenuation  than  shear 
waves  in  the  same  material. 
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2.4  Reflection  and  Transmission  of  Waves 

When  an  elastic  wave  hits  a  boundary  between  two  different  materials,  part  of  the  wave 
energy  is  reflected  and  part  of  it  is  transmitted  very  much  like  the  reflection  and 
transmission  of  light.  For  shear  waves,  the  force  and  the  tangential  displacement  have  to 
be  continuous  across  the  boundary,  and  for  longitudinal  waves,  the  particle  velocity  and 
the  acoustic  pressure  have  to  be  continuous.  Now  let’s  define  the  wave  impedance  Z  for 
the  media  at  both  sides  of  the  boundary: 

Z,  -c, [kg/m2s],  (14) 

with  /  being  the  index  indicating  the  material.  Both  p  and  c  are  material  dependent. 

The  impedance  is  different  for  shear  and  longitudinal  waves,  because  their  wave 
propagation  velocities  differ.  Most  likely,  the  impedance  is  complex,  because  the 
propagation  velocity  can  be  complex  (see  Equation  1 1).  When  the  material  in  which  the 
incident  wave  travels,  has  index  1  and  the  other  material  index  2,  then  from  the 
conditions  formulated  above  it  follows  that  the  reflection  coefficient  F  and  the 
transmission  coefficient  r  are 

A  7—7  A  27 

r  =  —  =  — — —  and  z  -  —  -  — — — ! —  with  1  +  r  =  r ,  (15) 

At  Z,  +  Z2  A ,  Z,  +  Z, 

with  A,  the  amplitude  of  the  incident  wave,  Ar  the  amplitude  of  the  reflected  wave,  and 
A,  the  amplitude  of  the  transmitted  wave  (see  Figure  3).  When  the  second  material  is 
very  dense  compared  to  the  first  material,  such  that  Z2  »  Z, ,  then  it  follows  from 
Equation  15  that  virtually  nothing  is  transmitted  and  that  the  wave  is  reflected 
completely.  The  amplitude  of  the  reflected  wave  is  the  negation  of  the  amplitude  of  the 
incident  wave.  In  other  words  a  phase  shift  of  n  rad  has  occurred.  When  Z,  »  Z2  full 
reflection  occurs  also  and  at  the  boundary  the  material  makes  large  movements. 

Another  special  case  is  when  Z2  =  Z, ,  for  instance  when  both  materials  are  the  same 
(not  necessarily).  As  expected,  the  incident  wave  is  fully  transmitted  and  nothing  is 
reflected  [Pain,  1993], 


A, 

A 

4 - : - 

medium  2 

A 

medium  1 

Z2 

A 

Figure  3  Reflection  and  transmission  of  a  wave  that  hits  a  boundary  between  medium  1  and  medium  2 
perpendicularly.  Z,  indicates  the  impedance  of  medium  /.  A,  is  the  amplitude  of  the  incident 
wave,  Ar  of  the  reflected  wave  and  A,  of  the  transmitted  wave. 

There  are  two  interesting  applications  when  we  have  two  materials  with  different 
impedances.  In  the  first  application,  we  want  to  get  optimal  transmission  of  waves  with 
as  little  reflection  as  possible.  In  the  second  application,  we  want  to  isolate  the  two 
materials  from  one  another.  In  other  words,  transmission  should  be  minimised.  Starting 
with  the  first  application,  the  impedances  of  both  materials  have  to  be  matched  by 
insertion  of  a  third  material  of  thickness  l  and  impedance  Zmawh  in  between.  Total 
transmission  occurs  when 
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yjz^Z-2  and  / 


^  match  ZltC malcfj 


CO 


(16) 


with  Aumldl  the  wavelength  of  the  wave  in  the  inserted  material.  One  can  see  that  the 

matching  impedance  and  the  thickness  are  different  for  shear  and  longitudinal  waves. 
Complete  transmission  only  occurs  for  one  wave  frequency  [Pain,  1993]. 

Isolation  requires  the  opposite,  namely  an  impedance  mismatch.  To  reduce  transmission 
a  material  has  to  be  added  in  between  that  has  an  impedance  that  differs  greatly  from 
the  impedances  of  the  materials  at  both  sides.  This  causes  almost  total  reflection.  When 
a  material  with  a  large  impedance  is  inserted,  it  needs  to  be  a  stiff,  heavy  material  and 
that  is  not  practical.  A  material  with  a  small  impedance  needs  to  have  a  resonance 
frequency  that  is  lower  than  the  wave  frequency  and  has  a  low  shear  modulus.  Care 
should  be  taken  that  there  are  not  unforeseen  rigid  connections  between  the  two 
materials  on  both  sides  of  the  inserted  layer  [Yerges,  1969]. 

The  concept  of  wave  reflection  and  transmission  can  be  extended  from  one  dimension 
to  two  dimensions.  Now  the  wave  number  as  well  as  the  propagation  direction  is  a 
vector.  When  a  wave  travels  in  a  positive  x-  and  y-direction  with  an  angle  9  from  the  x- 
axis,  the  solution  to  the  wave  equation  in  Equation  8  becomes 

u(x,y,t)  =  Ae1^-1'}  =  AeJ(<" -*«*(*)-#*  “(*» .  (17) 

When  this  wave  is  travelling  in  medium  1  and  it  strikes  the  boundary  with  medium  2, 
that  extends  itself  in  theyz-plane,  a  part  of  the  wave  is  reflected  under  an  angle  9r .  The 
relation  between  the  incident  and  reflected  angle  is  governed  by  Snell’s  law  of 
reflection 

er  =  ei,  (is) 


with  6 ,  the  angle  of  incidence.  Another  part  of  the  wave  is  transmitted  into  medium  2 
under  an  angle  9t ,  which  is  governed  by  Snell’s  law  of  refraction 


sin  9t  _  c,  _  k2 


sin  9t 


sin  9,  =  —sin  9  ■ 
k , 


(19) 


Both  Equation  18  and  19  are  solutions  of  the  well-known  Snell’s  relations 


sin  9t  _  sin  9r  _  sin  6 [ 


(20) 


The  refraction  and  transmission  of  an  incident  wave  are  depicted  in  Figure  4. 


Figure  4 


Reflection  and  transmission  of  a  wave  that  hits  a  boundary  between  medium  1  and  medium  2 
under  an  angle  0,  and  with  an  amplitude  A,.  The  wave  reflects  partly  under  an  angle  6r  and  with 
an  amplitude  Ar  and  it  is  partly  transmitted  under  an  angle  9,  and  with  an  amplitude  A,. 
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From  Equation  19  it  follows  that 


1  (kA 

M  -  sin2  0,  =±  , 

1-  — 

UJ 

sin2  #, 


(21) 


Equation  19  and  21  can  be  filled  in  Equation  17  to  obtain  the  wave  solution  in  medium 
2.  Only  the  positive  solution  of  Equation  21  results  in  a  physical  meaningful  wave 
solution.  The  amplitudes  of  the  reflected  and  refracted  waves  are  given  by  the  reflection 
coefficient  T  and  the  refraction  coefficient  r  : 


^  Ar  c-.cos6.~c.cos9, 

T  =  — -  -  — — — - - — — and 

A-  c2  cos  6,  +  c,  cos6l 


2 c2  cos  6, 


A,  c2  cos  6 \  +  c,  cos  6, 


with  1  +  T  =  r. 


(22) 


2.5  Surface  Waves 

An  interesting  phenomenon  occurs  when  in  Equation  21  >  k2  ;  the  wave  is  totally 

reflected  at  the  boundary,  sin#,  >  1 ,  and  cos#,  becomes  imaginary: 


\(k  V 

cos#,  =  ±jj  -i 
\K2j 


sin2  #,  - 1 . 


(23) 


When  Equation  19  and  23  are  substituted  in  Equation  17,  a  new  wave  solution  is 
obtained  that  equals 


u(x,y,t)  =  Ae 


cot  A  jxk2 


sin2  0,-\-yk2-^-s\x\0i 


(24) 


Only  the  negative  solution  of  Equation  23  is  valid.  From  Equation  24  it  becomes  clear 
that  a  different  kind  of  wave  than  mentioned  before  occurs.  It  travels  in  the  positive  y- 
direction  along  the  boundary.  The  amplitude  of  this  wave  attenuates  exponentially  in 
the  positive  x-direction,  creating  an  evanescent  field  in  medium  2.  This  kind  of  wave  is 
called  a  surface  wave  or  evanescent  wave  and  has  a  wave  number  ks  that  is  smaller 

than  k{  and  larger  than  k2 .  Therefore,  the  wave  propagation  velocity  of  the  surface 
wave  is  larger  than  the  velocity  in  medium  1 ,  it  is  smaller  than  the  velocity  in  medium 
2,  and  it  depends  on  9i  [Pain,  1993;  Cheng,  1989;  Petrashen,  2003], 

The  theory  for  surface  waves  above  does  not  take  into  account  the  difference  between 
the  wave  numbers  of  shear  waves  and  compression  waves.  To  get  a  more  complete 
picture  the  theory  can  be  extended:  When  either  a  compression  wave  or  a  shear  wave  is 
reflected  or  transmitted,  then  a  part  of  the  wave  energy  will  be  converted  into  a  shear  or 
compression  wave  respectively.  Snell’s  law  for  an  incident  longitudinal  wave  now 
becomes 
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sin#,,  sin#,.,  sin#,.v  sin#,,  sin#,  5 

C\  .1  C\.l  C\  ,.v  C2.l  C2..v 


(25) 


The  second  index  in  #v  indicates  whether  the  angle  belongs  to  a  shear  ( s )  or 

longitudinal  (/)  wave.  From  the  different  propagation  velocities  it  follows  that 
#, ,  =  Or ,  >  6r  5  and  #, ,  >  #,  t .  For  an  incident  shear  wave,  Snell’s  law  looks  very 

similar: 

sin#,  s  sin #r/  sin#,.v  sin#,,  sin#,  v 

— - = - - - —  = - —  =  — - •  (26) 

Cl,v  C\J  C\,s  C2 .1  C2  ,.v 

From  this  equation,  it  follows  that  #,  t  =  9r  s  <6rl  and  #, ,  >  #,  ( .  Shear  and 

longitudinal  surface  waves  occur  under  the  conditions  listed  in  Table  1 .  From  the  wave 
velocities,  it  follows  that  a  shear  surface  wave  always  is  accompanied  by  a  compression 
surface  wave.  In  addition,  a  shear  surface  wave  is  not  likely  to  occur  when  the  incidence 
wave  is  a  longitudinal  wave.  Most  likely  to  happen  is  a  longitudinal  surface  wave  when 
a  shear  wave  was  the  incident  wave  [Petrashen,  2003]. 


Table  1  The  conditions  for  which  a  longitudinal  or  shear  incident  wave  results  in  a  longitudinal  or  shear 
surface  wave.  It  is  always  true  that  ctj  >  cu  and  Cj.i  >  Ci.s  Therefore,  when  a  shear  surface  wave 
exists,  then  a  longitudinal  surface  wave  also  exists. 


surface  wave  occurrence 

longitudinal  shear 

incident  wave 

longitudinal 

sin  #, ,  >  — 

C2  J 

Ci  / 

sin  #, ,  >  — 1- 

’  C2,v 

shear 

sin  #,  v  >  —  or  sin  #,  v  >  — 

C2 ,/  c\ ,/ 

sin  #,  v  >  — 

"  C2,s 

2.6  Wave  Diffraction 

When  a  small  object  consisting  of  medium  2  is  embedded  in  medium  1,  waves 
propagating  in  medium  1  will  tend  to  bend  around  the  object.  This  is  called  diffraction. 
Diffraction  occurs  when  the  size  of  the  object  is  close  to  the  wavelength  of  the 
travelling  wave.  When  the  object  is  smaller,  then  the  wave  will  not  be  affected.  When 
the  object  is  larger,  then  reflection  and  transmission  will  occur. 

2.7  Finit  Sized  Vibrators 

Imagine  that  the  medium  is  excited  with  a  finite  sized  oscillating  plate  at  the  upper 
boundary.  Due  to  its  size,  the  waves  created  will  be  neither  pure  plane  waves,  nor 
spherical  waves,  but  somewhere  in  between.  As  a  result  the  waves  will  attenuate  more 
quickly. 

When  the  plate  oscillates  perpendicular  to  the  surface  of  the  medium,  compression 
waves  will  travel  downwards  and  shear  waves  along  the  surface.  However,  when  the 
plate  oscillates  over  the  surface  shear  waves  will  travel  downwards  and  compression 
waves  along  the  surface  (see  Figure  5).  This  will  occur  especially  when  the  oscillating 
plate  is  slightly  embedded  in  the  medium. 
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a  b 


Figure  5  a  When  the  vibrator  oscillates  horizontally,  a  shear  wave  moves  downwards  and  a 
longitudinal  wave  propagates  along  the  surface, 
b  When  the  vibrator  oscillates  vertically,  a  longitudinal  wave  moves  downwards  and  a  shear 
wave  propagates  along  the  surface.  A  combination  of  shear  and  longitudinal  waves  will 
occur  at  different  angles. 


2.8  Slip 

When  the  oscillating  plate  vibrates  along  the  surface,  shear  force  is  applied  to  the 
medium.  When  the  plate  is  not  rigidly  attached  to  the  underlying  medium  slip  can  occur 
and  wave  energy  is  lost.  How  easily  two  materials  slip  is  described  by  the  static  friction 
coefficient  fis ,  which  depends  on  the  normal  force.  If  the  shear  force  exceeds  the  static 
friction  coefficient,  the  material  will  slip  and  the  friction  force  will  reduce.  The  friction 
force  is  then  governed  by  the  kinetic  friction  coefficient  /uk  and  the  normal  force,  until 

the  shear  velocity  becomes  zero  [Weisstein,  2004a], 
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3  Structure  of  the  Skin 


3.1  Composition 

The  skin  has  a  thickness  range  of  0.5  to  6  mm  and  it  is  often  described  as  a  layered 
medium  of  two  distinctive  layers,  the  epidermis  and  the  dermis.  In  the  epidermis  (75- 
150  pm  thickness),  new  epidermal  cells  are  continuously  formed  and  transported 
outwards.  On  their  way,  these  cells  flatten  and  keratinise;  filaments  of  keratin  and  links 
between  these  filaments  are  formed  and  eventually  the  cells  die,  leaving  a  tough, 
protective  layer.  The  epidermis  becomes  more  dense,  stiffer  and  drier  towards  the 
outside. 

The  dermis  (0.6-3  mm  thickness)  is  the  thicker  layer  of  the  two  and  it  consists  mainly  of 
the  fibrous  proteins  collagen,  elastin  and  reticulin.  Together  they  form  a  complex 
network  filled  with  fluid.  This  network  embeds  blood  vessels,  nerves,  lymphatics,  hair 
follicles,  sweat  glands  and  mechanoreceptors.  The  papillary  layer  of  the  dermis  lies  just 
under  the  wavy  formed  junction  between  the  epidermis  and  the  dermis  (the  stratum 
basale)  and  consists  of  open  networks  of  fine  collagen  fibres.  Ninety  percent  of  the 
dermis  consists  of  a  denser  network  of  coarse  collagen  fibres  called  the  reticular  layer. 
This  means  that  the  reticular  layer  is  stiffer  than  the  papillary  layer.  In  relaxed  skin,  the 
collagen  fibres  lay  in  a  wave  pattern.  When  the  skin  is  being  stretched,  at  first  the 
collagen  fibres  will  align  allowing  the  skin  to  be  easily  stretched.  One  by  one  each  fibre 
will  stretch,  causing  a  high  stiffness  [Manschot,  1985;  Bischoff,  Arruda,  &  Grosh, 

2000].  Under  the  dermis  lies  a  layer  of  subcutaneous  fat  that  rests  on  a  stiff  layer  of 
muscle  or  bone. 

3.2  Mechanoreceptors 

The  different  types  of  mechanoreceptors  in  the  skin  work  together,  in  order  to  transmit 
tactile  information  to  the  brain.  These  mechanoreceptors  have  different  morphologies, 
functions  and  locations  (see  Figure  6).  Mechanoreceptors  in  the  human  hairy  skin  show 
some  similar  morphological  features  [Iggo  &  Andres,  1982]: 

-  The  receptor  is  formed  from  the  terminal  of  an  afferent  that  has  its  cell  body  in,  or 
near,  a  dorsal  root  ganglion  or  cranial  nerve  ganglion. 

-  The  receptor  terminal  shows  morphological  specialisation:  accumulation  of 
mitochondria;  presence  of  granular  vesicles;  fine  filamentous  ground  substance 
under  the  axonal  membrane. 

-  A  myelin  sheath  is  absent  from  the  receptor  terminal. 

-  The  receptor  terminal  often  is  connected  to  nonneural  cells,  which  may  form  a 
complex  and  characteristic  structure  in  which  the  nonneural  cells  can  act  as  a 
transducer. 

Hair  follicles  are  surrounded  by  mechanoreceptors  just  below  the  duct  of  the  sebaceous 
gland  if  one  is  present.  There  is  a  difference  between  the  innervations  of  guard  hairs 
(hairs  with  large  diameter  hair  shafts)  and  down  hairs  (smaller  hairs,  sometimes 
minute).  Guard  hairs  are  innervated  by  a  collar  of  lanceolate  terminals.  Each  terminal  is 
lance  shaped  and  sandwiched  between  two  Schwann  cells.  Small  axonal  spines  protrude 
from  the  sharp  sides  of  the  lance  into  the  surrounding  tissue.  The  lances  are  aligned 
parallel  to  the  long  axis  of  the  hair.  Together  they  form  a  palisade  around  the  hair. 
Lanceolate  terminals  are  believed  to  be  rapidly  adapting  [Munger  &  Ide,  1988;  Iggo 
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et  al.,  1982],  Several  stem  axons  supply  one  hair  follicle,  but  one  stem  axon  also 
supplies  several  hair  follicles  [Iggo  et  al.,  1982],  causing  overlap. 

Spiralling  around  the  palisade  of  lanceolate  terminals  is  the  Ruffini  terminal.  The 
terminal  is  not  encapsulated  by  Schwann  cells  as  opposed  to  Ruffini  terminals  found  in 
the  skin  near  joints  [Mungeret  al.,  1988;  lggo  et  al.,  1982;  Hamann,  1995].  To  make  a 
distinction  between  the  two  types  of  Ruffini  terminals,  the  naked  terminals  are  called 
piio-Ruffini  terminals.  The  Ruffini  endings  are  believed  to  be  slowly  adapting  receptors 
[Iggo  et  al.,  1982;  Hamann,  1995]  and  are  sensitive  to  vibrations  ranging  from  15- 
400  Hz  [Bolanowski,  Jr.,  Gescheider,  Verrillo,  &  Checkosky,  1988;  Johansson, 
Landstrom,  &  Lundstrom,  1982].  Opposed  to  guard  hairs,  down  hairs  often  lack  pilo- 
Ruffini  terminals  [Munger  et  al.,  1988],  but  are  surrounded  by  more  sensitive  lanceolate 
terminals  [Hamann,  1995]. 

In  between  hairs  clusters  of  Merkel  complexes  exist  [Munger  et  al.,  1988;  Iggo  et  al., 

1 982],  They  are  located  at  the  base  of  the  epidermis  and  consist  of  50  to  70  Merkel 
complexes  in  the  form  of  a  circular  or  oval  shape  with  a  size  of  approximately 
1 00  x  300  /jm  .  The  cluster  is  supplied  by  only  one  axon  [Iggo  et  al.,  1 982].  The 
epidermis  overlying  the  cluster  is  thickened  and  this  dome  is  called  a  Haarscheibe.  A 
Merkel  complex  consists  of  a  disc  shaped  axon  terminal,  called  a  Merkel  disc  that  is  in 
close  contact  with  a  special  cell,  called  Merkel  cell.  The  Merkel  cell  has  desmosomes 
with  its  adjoining  keratinocytes.  It  is  speculated  that  the  Merkel  cell  communicates  with 
the  Merkel  disc  through  synaptic-like  junctions.  Merkel  complexes  are  categorised  as 
SAI  receptors  and  they  are  sensitive  to  perpendicular  displacement  of  the  skin  [Munger 
et  al.,  1988;  Iggo  et  al.,  1982;  Hamann,  1995],  They  are  sensitive  to  vibrations  ranging 
from  0.4-100  Hz  [Johansson  et  al.,  1982;  Bolanowski,  Jr.  et  al.,  1988]. 

Deeper  in  the  skin,  in  the  dermis,  lie  the  largest,  most  sensitive  and  the  most  extensively 
investigated  mechanoreceptors,  called  the  Pacinian  corpuscles.  Pacinian  corpuscles 
consist  of  a  myelinated  axon  with  a  bare  end.  This  end  is  enveloped  by  20  to  70  densely 
packed  lamellae  consisting  of  an  inner  and  an  outer  core  [Hamann,  1995;  Iggo  et  al., 

1 982],  The  lamellae  form  a  oval  shape  with  a  length  of  0.5-2  mm  and  a  diameter  of  0.7 
mm  [Iggo  et  al.,  1982;  Munger  et  al.,  1988].  The  bare  axon  end  has  axonal  spines  of  6 
nm  that  protrude  between  the  layers  of  the  inner  lamellae,  very  much  like  the  axonal 
spines  found  in  lanceolate  terminals.  The  Pacinian  corpuscles  are  rapidly  adapting  and 
are  sensitive  to  vibrations  from  approximately  80  to  300  Hz  [Hamann,  1995;  Iggo  et  al., 
1982;  Munger  et  al.,  1988;  Bolanowski,  Jr.  et  al.,  1988;  Johansson  et  al.,  1982].  It  is 
believed  that  slip  in  the  lamellae  is  partly  responsible  for  the  high-pass  properties  of  the 
mechanoreceptor  [Munger  et  al.,  1988;  Hamann,  1995;  Sherwood,  1997]. 

In  great  abundance  in  the  skin  are  the  free  nerve  endings.  These  are  small  nerve  fibres 
without  specialised  nerve  endings.  Some  of  them  are  capable  of  detecting  mechanical 
stimuli.  They  might  give  a  tickling  sensation  when  stimulated  [Hamann,  1995].  The 
well  known  Meissner  corpuscles  are  not  present  in  hairy  skin. 
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Figure  6  Schematic  representation  of  the  distribution  of  mechanoreceptors  in  the  hairy  skin.  Shown  are 
lanceolate  terminals.  pilo-Ruffini  terminals.  Merkel  complexes,  and  Pacinian  corpuscles. 

3.3  Mechanotransduction  of  Vibrations 

Mechanotransduction  describes  the  process  needed  for  the  conversion  of  mechanical 
stimuli  into  a  bursts  of  action  potentials.  It  can  be  viewed  as  a  three-stage  process 
[French,  1992]:  firstly,  the  stimulus  is  mechanically  coupled  to  the  receptor.  Some 
responsible  structures  already  mentioned  above  are  hairs,  lamellae,  axonal  spines,  and 
the  epidermis.  These  structures  attenuate  or  amplify  the  stimulus,  acting  like  a  specific 
filter.  Membrane  penetrating  proteins  may  have  a  function  in  mechanotransduction.  An 
important  family  of  proteins  could  be  degenerin/epithelial  Na+  channel  proteins  [Welsh, 
Price,  &  Xie,  2002], 

Secondly,  the  deformation  of  the  receptor  is  converted  into  an  electrical  signal  by 
mechanically  sensitive  ion  channels.  These  channels  appear  to  be  tension  sensitive 
[French,  1992],  Interesting  are  the  results  of  experiments  Hamill  performed  on 
mechanosensitive,  nonsensory  frog  oocytes  (eggs)  [Flamill  &  McBride,  1995].  Hamill 
applied  suction  to  the  cell  membranes  of  the  oocytes  and  measured  the  resulting  change 
in  membrane  current  with  a  patch  clamp.  When  a  pressure  ramp  was  applied  on  a  piece 
of  membrane  containing  only  one  channel,  the  membrane  current  showed  a  step 
response.  When  more  channels  were  involved  the  current  increased  with  increasing 
pressure  and  saturated  before  the  maximum  pressure  was  reached.  Flamill  concluded 
from  this  that  channels  are  either  fully  opened  or  closed  and  that  each  channel  opens  at 
a  different  pressure.  He  also  applied  a  step  stimulus  and  measured  the  adaptation  of 
multiple  channels  and  single  channels.  It  turned  out  that  each  channel  had  a  different 
adaptation  time,  resulting  in  exponential  adaptation  for  a  group  of  channels. 
Remarkably,  adapted  channels  could  be  reactivated  by  further  increasing  the  pressure. 
Thirdly,  the  graded  potential  (depolarisation)  caused  by  the  influx  and  efflux  of  ions  has 
to  lead  to  action  potentials.  When  the  membrane  potential  in  the  axon  near  the  receptor 
reaches  a  threshold  potential,  voltage  controlled  ion  gates  open  and  a  sudden  influx  of 
Na+  occurs  resulting  in  a  sharp  increase  of  the  local  membrane  potential 
(depolarisation).  The  Na+  influx  is  rapidly  followed  by  an  efflux  of  K+,  resulting  in  a 
sharp  decrease  of  the  membrane  potential  (repolarisation).  Together  they  form  the 
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potential  peak  characteristic  for  the  action  potential  [Sherwood,  1997],  Several 
mechanisms  can  influence  the  graded  potential  to  action  potential  encoding,  which  are 
described  by  French  [Welsh  et  al.,  2002]. 

3.4  Wave  Propagation  Models  for  the  Skin 

In  the  literature,  several  mechanical  models  of  skin  have  been  described.  Usually  they 
describe  the  skin  as  being  composed  of  one  to  three  layers  of  materials  with  different 
mechanical  properties.  Below  some  models  and  results  are  described.  Mechanical 
parameters  are  listed  when  relevant  or  present. 

Sherrick  and  Cholewiak  describe  the  skin  as  a  ‘semi  plastic  membrane  that  overlays  a 
viscoelastic  substrate  that  varies  greatly  in  local  density,  hardness,  and  other  physical 
properties’.  On  the  skin  surface,  mechanical  disturbances  propagate  as  a  continually 
damped  wave,  which  attenuates  following  the  inverse-square  law.  The  surface  wave  is 
generated  by  a  combination  of  shear  and  compression  waves  within  the  substrate 
[Sherrick  &  Cholewiak,  1986]. 

Von  Gierke  et  al.  recognised  three  types  of  waves  that  propagate  through  the  skin: 
compression  waves,  shear  waves  and  surface  waves.  According  to  Von  Gierke  shear 
waves  attenuate  much  faster  in  the  skin  than  compression  waves.  As  a  model,  he  used  a 
vibrating  rigid  sphere  in  an  unlimited,  homogeneous,  viscous,  elastic,  compressible 
medium.  The  parameter  values  he  has  used  are  listed  in  Table  2.  Note  that  the 
imaginary  parts  of  the  Lame  constants  are  linearly  dependent  on  the  frequency  [von 
Gierke,  Oestereicher,  Franke,  Parrack,  &  von  Wittern,  1952], 

Table  2  Mechanical  parameter  values  as  used  by  Von  Gierke  [von  Gierke  et  al.,  1952], 


p 

P  =  P  +  JUP 

A  —  X  +  j  coA 

M 

P 

X 

X 

kg/m3 

Pa 

Pas 

Pa 

Pas 

1.1  to3 

2.5  103 

15 

2.6- 109 

0 

Ashby  et  al.  list  some  mechanical  properties  for  a  large  number  of  natural  materials. 
These  natural  materials  also  include  skin.  For  the  skin  they  had  taken 
p  =  1 .3  •  1 03  kg/m3  and  E  =  40  •  1 06  Pa  [Ashby,  Gibson,  Wegst,  &  Olive,  1 995]. 
Van  Doren  modelled  the  skin  of  the  fingertip  as  a  slab  of  homogeneous,  isotropic, 
viscoelastic  material  with  the  material  properties  listed  in  Table  3.  The  slap  was  rigidly 
connected  to  the  ground.  He  measured  the  spatiotemporal  sensitivity  at  two  temporal 
frequencies  (8  Hz  and  128  Hz)  and  twelve  spatial  frequencies  ranging  from  0.0-1 .03 
cycles/mm.  The  results  showed  a  decrease  in  the  sensitivity  with  increasing  spatial 
frequency  for  128  Hz,  and  an  increase  in  sensitivity  for  8  Hz.  He  assumed  that  the 
psychophysical  threshold  of  a  receptor  type  corresponds  to  a  constant  value  of  one 
strain  component  evaluated  at  the  location  of  the  appropriate  receptor.  He  used  the 
model  of  the  skin  and  the  measurements  to  obtain  this  strain  component  and  the 
location  of  the  mechanoreceptor  in  the  skin.  The  measurements  were  also  used  to 
estimate  the  skin  thickness  in  the  model  [Doren,  1989]. 
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Table  3  Mechanical  parameter  values  as  used  by  Van  Doren  [Doren,  1989], 


Medium 

Layer 

P 

P  =  P  +  jcop' 

X 

index 

thickness 

p 

P 

m 

kq/m3 

Pa 

Pa.s 

Pa 

1 

4.96-10'3 

1.1  to3 

2.5-104 

40 

oo 

Klochkov  modelled  the  skin  as  a  two-layer  medium.  The  first  layer  has  a  thickness  that 
is  varied.  The  second  layer  extends  into  infinity.  He  had  positioned  a  surface  source  of 
forced  vibration  on  top  of  the  first  layer,  capable  of  normal  and  tangential  oscillations. 
For  normal  displacement,  the  frequency  dependency  was  calculated.  It  appeared  that  the 
resonance  frequency  increased  with  a  decrease  in  the  layer  thickness  and  that  the 
resonance  peak  became  lower  and  wider.  The  resonance  frequencies  ranged  from  50  to 
300  Hz.  In  the  case  of  tangential  displacement,  two  resonance  frequencies  exist  with  a 
dip  in  between.  For  layer  thicknesses  larger  than  2  cm  the  amplitude  of  the  dip  became 
less  and  its  frequency  shifted  to  lower  values.  For  layer  thicknesses  smaller  than  2  cm 
the  frequency  of  the  dip  shifted  to  higher  values  and  the  amplitude  decreased.  In  the 
case  of  a  layer  thickness  of  2  cm  the  two  resonance  frequencies  were  approximately  45 
and  90  Hz.  The  dip  occurred  at  56  Hz.  Klochov  also  calculated  the  two  dimensional 
(along  the  surface  and  in  the  depth  of  the  tissue)  spatial  distribution  of  the  amplitude  of 
the  displacement  for  both  tangential  and  normal  oscillations.  The  oscillation  frequency 
was  set  to  70  Hz.  The  results  showed  that  in  the  case  of  normal  oscillation  the  waves 
damped  out  quickly  in  both  dimensions:  within  5  cm  along  the  surface,  and  within  3  cm 
in  the  depth.  Tangential  oscillation  resulted  in  a  more  complex  distribution  that  showed 
less  attenuation  [Klochkov,  2002). 


Table  4  Mechanical  parameter  values  as  used  by  Klochkov  [Klochkov,  2002]. 


Medium 

index 

Layer  thickness 

p 

P  =  P  +  jap" 

X 

P 

P 

m 

kq/m3 

Pa 

Pa.s 

Pa 

1 

0.3-4-10'2 

1.05-103 

5-103 

2.5 

2.67- 109 

2 

OO 

1.6-103 

6.4- 109 

0 

5.7- 109 

Pereira  modelled  the  skin  as  a  material  with  three  layers:  the  stratum  comeum,  the 
dermis,  and  the  subcutaneous  fat.  On  top  of  the  stratum  comeum,  he  placed  a 
tangentially  oscillating  line  source.  He  varied  the  shear  modulus  and  the  loss  modulus, 
while  measuring  the  surface  wave  velocity.  Note  that  the  loss  modulus  is  independent 
of  the  frequency.  Pereira  concluded  that  the  stiffness  of  the  stratum  comeum  had  a  large 
influence  on  the  propagation  velocity  of  surface  waves  [Pereira,  Mansour,  &  Davis, 
1991). 
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Table  5  Mechanical  parameter  values  as  used  by  Pereira  [Pereira  et  al..  1991],  Medium  1  is  the  stratum 
comeum.  medium  2  is  the  dermis,  and  medium  3  is  subcutaneous  fat. 


Medium 

Layer  thickness 

P 

P  =  P  +  jp" 

index 

P 

P 

m 

kg/m3 

Pa 

Pa 

1 

15-10-6 

1-103 

IOTOV2 

0.25-0.75.^7 

2 

15  10  3 

1-103 

60- 103 

32.5-103 

3 

OO 

0.9-103 

20- 103 

10.103 

Timanin  also  used  a  three-layered  model  with  a  vibrating  pressure  source  on  top.  He 
used  excitation  frequencies  ranging  from  20  to  60  Hz.  Lower  frequencies  resulted  in 
less  attenuation.  The  displacement  amplitude  at  a  fixed  point  increased  with  the  radius 
of  the  pressure  source  [Timanin,  2002]. 


Table  6  Mechanical  parameter  values  as  used  by  Timanin  [Timanin,  2002],  The  2’s  have  been 
estimated  using  the  given  longitudinal  wave  propagation  velocities  c/. 


TNO  report  |  TNO-DV3  2005  A075 


23/29 


4  Design  Criteria 


The  knowledge  from  the  previous  two  chapters  can  be  used  to  design  a  tactile  display 
that  is  capable  of  precise,  local  skin  stimulation.  In  this  chapter  some  design  criteria  are 
composed.  For  this  purpose  the  tactile  display  has  been  divided  into  three  parts:  the 
tactor,  the  tactor-skin  interface,  and  the  part  that  holds  the  tactors,  which  has  been 
labelled  as  the  tactor  carrier.  For  each  of  these  parts  criteria  listed  below. 

4.1  Tactor 

The  correct  answer  to  the  question:  ‘Which  type  of  mechanoreceptor  is  stimulated  by 
mechanical  stimulation  on  the  surface  of  the  skin?’  is:  all  types.  Flow  it  is  perceived 
depends  on  the  location,  sensitivity,  density  and  the  type  of  innervation  of  a  type  of 
mechanoreceptor.  Mechanoreceptors  located  near  the  surface  are  subjected  to  a  larger 
wave  amplitude  then  the  more  deeply  located  mechanoreceptors.  Considering  this,  the 
Pacinian  corpuscles  have  a  disadvantage  here,  because  they  are  located  deeply  in  the 
skin.  In  addition  the  density  of  Pacinian  corpuscles  in  the  skin  is  very  low.  However, 
they  are  believed  to  be  very  sensitive,  because  they  are  large  and  have  many  layers  that 
effectively  transfer  vibration  energy  to  the  axonal  spikes.  Mechanoreceptors  that  are 
located  near  the  surface  and  that  are  densely  packed  in  the  skin  are  the  clusters  of 
Merkel  complexes  in  between  hairs.  The  localisation  of  a  stimulus  might  not  be  very 
good,  because  each  cluster  is  supplied  by  only  one  axon,  whereas  one  Pacinian 
corpuscle  is  connected  to  a  single  axon.  A  more  promising  innervation  have  the 
lanceolate  terminals  in  hair  follicles;  one  hair  follicle  is  supplied  by  more  than  one 
nerve  fibre.  Pilo-Ruffini  endings  probably  also  play  a  mayor  role  in  the  detection  of 
vibration  stimuli  on  locations  where  guard  hairs  grow.  It  is  expected  that 
mechanoreceptors  are  best  stimulated  by  shear  waves,  because  according  to  Equation  9 
shear  waves  exert  a  local  force  difference  on  two  neighbouring  particles  that  is 
approximately  a  thousand  times  larger.  Consequently,  shear  stimulation  of  the  skin 
results  in  more  local  stress  in  the  mechanoreceptors;  the  receptor  is  stretched  and 
compressed  instead  of  pushed  forwards  and  backwards  when  a  wave  passes.  However, 
a  vibrator  that  oscillates  along  the  surface  of  the  skin  slips  more  easily  than  a  vibrator 
that  oscillates  perpendicular.  This  results  in  a  less  efficient  transfer  of  the  vibration 
energy.  Vibrators  that  are  used  in  a  tactile  display  often  do  not  have  a  clear  oscillation 
direction  and  will  evoke  a  mixture  of  longitudinal  and  shear  waves  with  varying 
polarisations. 

Since  the  epidermis  is  stiffer  than  the  dermis,  the  waves  will  propagate  farther  in  the 
epidennis  parallel  to  the  surface  of  the  skin.  According  to  Equation  15,  shear  and 
longitudinal  waves  will  most  likely  penetrate  through  the  dermis  and  the  subcutaneous 
fat  layer  and  reflect  almost  entirely  on  the  much  stiffer,  underlying  layer  that  consists  of 
muscle  or  bone.  The  reflected  waves  arrive  at  the  epidermis-air  boundary  under  every 
possible  angle,  because  of  the  spherical  like  form  of  the  waves.  The  propagation 
velocity  of  longitudinal  waves  in  air  is  smaller  than  the  velocity  in  skin  and  shear  waves 
cannot  exist  in  air.  Therefore,  according  to  Table  1  only  longitudinal  surface  waves  that 
originate  from  shear  waves  can  exist.  Perpendicular  excitation  of  the  skin  would  of 
course  result  in  a  shear  wave  travelling  along  the  surface  of  the  skin  starting  at  the  rim 
of  the  source,  because  of  the  up  and  down  motion  of  the  skin.  In  the  case  of  parallel 
excitation,  a  longitudinal  wave  would  travel  along  the  surface  (see  Figure  5). 

The  attenuation  of  waves  in  the  skin  depends  on  the  loss  modulus  //  and  is  governed 
by  Equation  13.  This  attenuation  increases  with  the  vibration  frequency.  From 
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Table  2-6  it  becomes  clear  that  A  for  the  dermis  is  a  factor  106  larger  than  ju  ,  which 
results  in  a  much  smaller  attenuation  for  longitudinal  waves  than  for  shear  waves.  It  can 
be  calculated  that  a  100  Hz  plane  shear  wave  has  to  travel  approximately  4  cm  in  the 
dermis  to  become  a  factor  e  smaller.  Given  that  the  thickness  of  the  skin  is  often  less 
then  6  mm  and  that  in  general  tactors  are  relatively  large  compared  to  the  thickness  of 
the  skin,  it  can  be  concluded  that  directly  underneath  a  tactor  at  a  wave  frequency  of 
100  Hz  mechanoreceptors  experience  almost  the  same  wave  amplitude. 

Diffraction  could  cause  high  attenuation  of  waves,  because  of  increased  path  lengths 
and  complicated  interference  patterns.  However,  the  propagation  velocity  of  shear 
waves  in  the  dermis  will  be  approximately  cs  »  2  m/s  and  the  wavelength  at  100  Hz 

will  be  approximately  2  cm  (for  longitudinal  waves  they  are  both  approximately  a 
thousand  times  larger).  Structures  in  the  skin  like  cells,  nerves,  vessels,  hair  follicles 
and  glands  are  much  smaller  than  that.  Therefore,  shear  waves  and  certainly 
longitudinal  waves  will  not  diffract  in  the  skin.  When  the  wave  frequency  is  high 
enough  (/  »  100  Hz  )  scattering  will  occur.  However,  these  frequencies  lie  above  the 
frequencies  that  can  be  perceived  by  the  mechanoreceptors  in  the  human  skin. 

Waves  with  wave  frequencies  larger  than  100  Hz  will  not  stimulate  the  slowly  adapting 
Merkel  cells,  but  they  will  stimulate  the  very  sensitive  Pacinian  corpuscles. 
Furthermore,  they  dampen  out  more  quickly  along  the  surface  of  the  skin  then  slowly 
oscillating  waves.  Although  vibrators  often  will  use  more  energy  when  operating  at 
higher  frequencies,  it  is  preferable  to  use  a  stimulation  frequency  of  200  Hz  (within  the 
frequency  range  in  which  Pacinian  corpuscles  are  the  most  sensitive). 

Table  7  Summary  of  the  design  criteria  for  the  tactor. 

Design  criteria  for  the  tactor _ _ _ _ _ _ _ 

When  no  slip  across  the  skin  can  occur,  use  a  tactor  that  oscillates  parallel  to  the  skin  for  better 

stimulation  of  the  mechanoreceptors.  Otherwise  use  perpendicular  stimulation. 

Use  a  stimulation  frequency  of  200  Hz 


4.2  Tactor-Skin  Interface 

The  interface  between  the  tactor  and  the  skin  should  be  designed  such  that  the  transfer 
of  vibration  energy  from  the  tactor  to  the  skin  is  large.  A  well  designed  contact  will 
reduce  the  need  for  strong  vibrators  that  have  high  energy  consumption. 

Equation  16  tells  us  to  insert  a  material  between  the  vibrator  and  the  skin  with  an 

impedance  that  equals  -JZ^Z^  and  a  thickness  of  A/ 4  .  A  perfect  match  is  not  feasible 

in  the  case  of  longitudinal  waves,  because  its  wavelength  in  solid  materials  is  meters 
long.  For  shear  waves,  the  matching  of  impedances  will  be  easier.  However,  the 
impedance  of  the  skin  varies  with  the  location  on  the  skin  so  the  match  will  never  be 
perfect. 

Another  way  to  improve  the  transfer  of  wave  energy  is  achieved  by  gluing  the  tactor  to 
the  skin.  This  will  give  the  highest  friction  between  the  tactor  and  the  skin.  Often  it  will 
be  impractical  to  attach  the  vibrator  rigidly  to  the  skin.  In  that  case,  the  vibrator  will  be 
pressed  onto  the  skin  with  or  without  a  material  like  clothing  in-between.  In  the  first 
case,  the  vibrator  should  have  enough  static  friction  or  else  it  will  slip.  In  the  latter  case, 
it  is  expected  that  the  fibres  in  clothes  slip  easily  and  reduce  the  energy  transfer.  With 
simulation,  it  should  be  possible  to  predict  the  influence  of  different  types  of  clothing 
between  the  tactor  and  the  skin.  Mechanical  properties  of  common  textiles  are  listed  in 
Appendix  A. 
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The  friction  between  the  tactor  and  the  skin  can  also  be  increased  by  applying  more 
pressure  on  the  tactor.  However,  this  will  stretch  the  skin  and  since  it  is  a  nonlinear 
material,  it  becomes  stiffer  when  it  is  stretched  (larger  fj,  ).  This  is  caused  by  the 
mechanical  properties  of  collagen  networks.  In  a  stiffer  material,  there  will  be  less 
attenuation  of  waves  and  a  higher  wave  propagation  velocity.  Skin  stretch  can  also 
occur  during  body  movement  or  during  contact  with  objects  in  the  environment  or 
during  stimulations  with  large  amplitudes.  To  keep  a  constant  sensation,  the  amplitude 
of  pressure  fluctuation  have  to  be  reduced.  This  could  be  accomplished  by  letting  the 
force  act  on  a  large  surface. 

The  interface  between  the  tactor  and  the  skin  also  plays  a  role  in  localising  the 
vibrations  in  the  skin;  All  waves  will  attenuate  with  the  distance  from  the  source, 
obeying  approximately  Equation  10.  When  the  contact  area  between  the  vibration 
source  and  the  skin  is  large,  the  waves  will  behave  more  like  plane  waves.  When  it  is 
small,  the  wave  will  behave  more  like  spherical  waves.  Consequently,  vibrations 
originating  from  a  small  contact  area  will  dampen  more  quickly  than  waves  originating 
from  a  large  contact  area.  A  small  contact  area  is  preferable,  since  it  reduces  the 
spreading  of  the  vibrations.  It  also  reduces  the  need  for  strong  vibrators,  because  less 
skin  mass  has  to  be  moved. 

Table  8  Summary  of  the  design  criteria  for  the  tactor-skin  interface. 

Design  criteria  for  the  tactor-skin  interface _ 

Glue  the  tactor  to  the  skin,  or  coat  the  tactor  with  a  material  that  does  not  slip  easily  to  prevent 
slip _ 

Keep  the  tactor-skin  contact  area  small  in  order  to  keep  the  vibrations  local _ 

Insert  material  of  thickness  X/4  between  the  tactor  and  the  skin  (if  possible)  to  obtain  a  good 
transfer  of  vibration  energy _ 


4.3  Tactor  Carrier 

In  the  current  version  of  the  tactile  vest  columns  tactors  are  attached  to  strips  of  Velcro. 
When  one  tactor  in  a  column  is  activated,  waves  propagate  along  the  strip  of  stiff 
Velcro.  These  waves  transfer  energy  to  other  tactors  and  introduce  tactile  noise  in  the 
skin.  A  vibrator  can  be  isolated  from  its  carrier  by  embedding  it  in  a  soft  material  with 
an  impedance  that  greatly  differs  from  the  impedance  of  the  vibrator,  which  on  its  turn 
is  surrounded  by  another  dense  material.  This  second  layer  causes  the  vibration  to  be 
reflected  back  to  the  source.  Care  should  be  taken  that  vibrations  do  not  leak  to  the  bulk 
material  via  the  skin  (see  Figure  7).  Fixation  of  the  skin  just  around  the  vibrator  will 
damp  the  waves  travelling  on  and  near  the  surface  of  the  skin,  but  not  the  waves  deeper 
below  the  surface.  Electric  wires  can  also  create  vibration  leaks. 
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Figure  7  Proposal  for  embedding  of  the  vibrator.  The  transfer  of  vibration  energy  from  the  vibrator  to 
the  skin  is  maximized  by  the  addition  of  an  extra  layer  of  material,  called  the  matching 
material.  The  transfer  of  vibration  energy  from  the  vibrator  to  the  bulk  material  is  minimised  by 
a  layer  of  soft  material  followed  by  a  layer  of  dense  material. 

In  some  user  environments,  a  low  noise  emission  is  important.  The  amount  of  noise 
produced  by  a  vibrator  depends  on  the  transfer  of  vibration  energy  to  the  air.  The 
isolation  of  the  vibrator  from  the  bulk  material  mentioned  above  reduces  this  transfer. 

In  general,  the  size  of  the  vibrating  surfaces  should  be  kept  as  small  as  possible. 

Table  9  Summary  of  the  design  criteria  for  the  tactor  carrier. 

Design  criteria  for  the  tactor  carrier _ _ _ _ _____ 

Keep  the  size  of  vibrating  surfaces  as  small  to  reduce  noise 

Pack  the  vibrator  in  isolating  material  to  prevent  spreading  of  the  vibrations  through  the  tactor 
carrier 

Do  not  attach  the  tactors  directly  to  dense  materials  in  the  tactor  carrier 
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A  Mechanical  Properties  of  Materials 


In  this  appendix  mechanical  properties  of  the  skin  and  common  clothing  are  given,  in 
case  future  tactile  systems  are  going  to  be  modelled  with  for  instance  finite  element 
(FEM)  models. 

Unfortunately,  it  was  not  possible  to  find  mechanical  properties  of  viscoelastic  skin, 
therefore  estimations  have  been  made  from  Table  2-6. 


Table  A.  1  Ranges  of  the  mechanical  properties  of  human  skin  estimated  from  Table  2-6. 


Material 

Density 

Shear 

Shear 

Volume  compres- 

Volume 

name 

elasticity 

viscosity 

sibility 

viscosity 

P 

M’ 

M" 

K 

K' 

kg/m3 

Pa 

Pas 

Pa 

Pa  s 

•103 

■103 

•109 

epidermis 

0.9-1. 0 

1.0 

6 

2.0 

0 

dermis 

1.0-1. 1 

2.5-5 

2.5-40 

2.6 

0 

Subcuta¬ 
neous  fat 

1.1 

1.7 

2 

1.0 

0 

In  the  textile  industry  the  Kawabata  test  method  is  used  to  obtain  mechanical 
parameters.  Unfortunately,  these  parameters  are  not  clearly  related  to  the  parameters 
used  in  this  document.  Important  Kawabata  parameters  are  the  fabric  weight  per  unit 
area  (fV),  the  fabric  thickness  ( T0 )  and  the  tensile  energy  ( WT  ).  WT  is  defined  as  the 
area  under  the  stress-strain  curve  from  0  to  maximum  strain  (maximum  stress  equals 
4.9N/cm  ).  Also  important  is  the  linearity  of  load-extension  curve  ( LT ),  which  is 
defined  as 


LT  = 


WT 

■j(max(f)-  max(cr)) 


(Al) 


Other  parameters  are  the  compressional  energy  ( WC )  and  the  linearity  of  the 
compression  thickness  curve  (  LC ).  They  are  similar  to  WT  and  LT  .  The  maximum 
stress  in  the  stress-strain  curve  equals  0.49  N/ cm2  .  Another  important  parameter  is 
the  shear  rigidity  ( G )  along  the  surface  obtained  from  the  shear  stress  versus  the  shear 
angle  curve.  This  curve  is  measured  by  shearing  a  textile  sample  of  dimensions  5  by  20 
cm  parallel  to  its  long  axis,  while  keeping  a  constant  tension  of  98.1  mN/cm .  The 
parameter  2  HG  gives  the  hysteresis  width  of  this  curve  at  the  shear  angle  of  0.5 
degrees. 

Obviously,  W  is  related  to  p  .  It  might  be  possible  to  make  estimations  of  the  Young’s 

moduli  at  zero  strain  along  the  surface  of  the  textile  and  perpendicular  to  the  surface  of 
the  textile  from  the  parameters  WT ,  LT ,  WC ,  and  LC  : 


max((j)2 
2  ■  WT  ’ 


(A2) 


in  which  E\ 


is  the  estimated  Young’s  modulus  around  zero  strain.  LT  and  ITT  can 


n-o 

be  replaced  by  LC  and  WC  respectively. 
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The  parameter  G  might  be  used  to  estimate  the  shear  modulus  and  the  imaginary  part 
of  the  shear  modulus  might  be  (only  partially)  obtained  from  2 HG .  Ranges  of  values 
for  the  parameters  mentioned  above  are  listed  in  Table  A. 2. 


Table  A. 2  Ranges  of  values  of  the  Kawabata  parameters  for  common  textiles. 


Parameter  name 

Value 

W 

60-525 

g/m2 

To 

0.4-3. 1 

mm 

WT 

5-50 

J/m 

LT 

0.3-0.95 

WC 

0.15-2.5 

J/m2 

LC 

0.3-0. 7 

G 

0. 5-8.0 

N/m° 

2HG 

3-20 

N/m 
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15.  ABSTRACT  (MAXIMUM  200  WORDS  (1044  BYTE)) 

When  tactors  of  a  tactile  display  are  mounted  in  a  tactor  carrier,  e.g.  a  vest  or  a  chair,  vibrations  will  spread  out 
through  the  skin  and  the  carrier.  Vibrations  should  be  local  to  ensure  good  detection;  spreading  through  the  skin  or 
the  environment  has  to  be  minimised. 

A  literature  study  has  been  done  to  obtain  the  following  design  criteria:  Tactors  should  oscillate  parallel  to  the  skin  at 
a  frequency  of  200  Hz.  Tactors  should  be  glued  to  the  skin  or  coated  with  a  material  that  does  not  slip  easily.  When 
the  tactors  cannot  be  attached  free  of  slip,  use  tactors  that  oscillate  perpendicular  to  the  skin.  The  size  of  the  contact 
area  between  the  tactor  and  the  skin  should  be  small.  Soft  material  encased  in  stiff  material  should  be  used  between 
the  tactor  carrier  and  the  tactor  in  order  to  reduce  the  amount  of  vibration  energy  emitted  to  the  tactor  carrier  and  to 
reflect  energy  to  the  skin. 
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